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Abstract 

We show how the homogeneous variational bicomplex provides a useful formalism 
for describing a number of properties of single-integral variational problems, and 
we introduce a subsequence of one of the rows of the bicomplex which is locally 
exact with respect to the variational derivative. We are therefore able to recover a 
Lagrangian from a set of equations given as a variationally-closed differential form. 
As an example, we show how to recover a first-order Lagrangian from a suitable set 
of second-order equations. 
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1 Introduction 



I met Willy Sarlet in 1986 in Gent, at the first of the annual Workshops on 
Differential Geometric Methods in Mechanics; since then we have had many 
discussions on topics of mutual interest, and worked together (along with 
Frans Cantrijn) on a series of papers on non-holonomic mechanics. The present 
paper is on a rather different subject, one which I have been studying on and 
off for some time. In fact I did not give a talk at the meeting to honour 
Willy's sixtieth birthday, as he had already heard parts of this paper twice 
before and I did not feel that it would be appropriate to trouble him with it 
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again. But at the meeting itself there was some discussion of the equations 
to be satisfied when 'Euler-Lagrange forms' arise from a Lagrangian in the 
fibred case, and so I was prompted to see how these equations might work in 
the homogeneous context. These new ideas are based on the constructions I 
had described previously, and can be expressed in a manner which is perhaps 
more elegant than in the fibred case. I shall restrict attention here to the case 
of a single-integral variational problem, although much of the discussion may 
be extended, albeit in a significantly more complicated manner, to multiple- 
integral problems. 

My interest in homogeneous variational problems was encouraged by Mike 
Crampin, and some aspects of this work, particularly that described in Sec- 
tion^ may be found in three joint papers published a few years ago [2ll3lPf] ; 
these papers also cover multiple-integral problems. 



2 Homogeneous variational problems 

There is an argument for saying that variational problems fall into two distinct 
types, depending upon the nature of the solutions. In one type, the solutions 
will be parametrised submanifolds of the ambient manifold; in the other, they 
will be submanifolds without any preferred parametrisation. 

Different geometrical structures are used to describe these two types of prob- 
lem. In the first case, we would typically have a fibration tt : E — > M and 
consider a Lagrangian on a jet manifold J k Tc; the solutions are local sections 
<p of E — > M, and the solution submanifolds lm(0) are parametrised by the 
sections themselves. Problems of this type arise in mathematical physics: for 
instance, in classical mechanics, where typically M = R, and in field theories, 
where dimM > 1. 

The other type of problem is that of minimal submanifolds, and in the first- 
order single-integral case this is just Finsler geometry. Here the ambient man- 
ifold E has no fibration, and instead of the affine jet manifolds J k n we could 
use instead the manifolds of m-dimensional contact elements J k (E,m) (oth- 
erwise known as jets of immersions, or as higher-order Grassmannians) ; the 
particular case J l {E, 1) is just the projective tangent bundle of E. But another 
approach, described in our papers [2l3fl] and elsewhere, is to use instead the 
manifolds of non-degenerate m-velocities T k m E. These are the total spaces of 
principal bundles T k m E — > J k (E, m) with structure groups L k m , the fc-th order 
m-dimensional jet groups; using these manifolds means that we should con- 
sider homogeneous objects, with a specific behaviour under the group action. 
There is an argument that working directly on J k (E,m) is more straightfor- 
ward, as there is no need to check invariance. But actually that is a disadvan- 
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tage, as starting with homogeneous objects one often finds that the resulting 
constructions are not invariant under the group action, and indeed this fact 
explains some of the difficulties found in the case of higher-order field theories. 

There is, of course, a relationship between the two types of problem. The 
'homogenisation trick' has been known for many years in classical mechanics: 
here, one takes a Lagrangian L(t, q l , q l ) given as a function with respect to the 
standard volume form dt, and replaces it by a homogeneous Lagrangian 

L(t,i,q i ,q i )=iL(t,q i ,q i i- 1 ) 



where t is now regarded as a variable with the same status as the q % . Geomet- 
rically, one has moved from the jet manifold </% (where 7r : E — > M. is some 
fibred manifold) to a suitable open subset of the slit tangent manifold T Q E, 
essentially by taking the contraction of L dt with the total time derivative. 
Much the same sort of trick works for higher orders, and also for multiple 
integral problems. 

In this paper we shall concentrate on the case m = 1, and so consider the 
higher-order non-degenerate tangent bundles T^E. We take as known (see, for 
instance, [5]) the existence of two canonical constructions on these manifolds. 
The first of these is the total time derivative g?t, a vector field along the 
projection T k+1 E — > T^E expressed in coordinates as 

k ■ d 

rf T = Ei+i)oT- 

p=0 °%) 



where q 1 are local coordinates on E and qL\ are the corresponding derivative 
coordinates with 'p dots' on a tangent manifold. The second is the vertical 
endomorphism S, a type (1, 1) tensor field on T^E expressed in coordinates 



as 

fc-i 







s = £( P +i)d4&— — 



In order to express conditions of homogeneity, we need to make use of some 
more vector fields A p on T^E; these are defined as the contractions A p = 
S p (dT) for 1 < p < k. Note that, although g?t is a vector field along a projec- 
tion, the properties of S ensure that A p are genuine vector fields on a man- 
ifold: they are, indeed, the fundamental vector fields of the principal bundle 
T k E _^ J k (E,l). 

A Lagrangian in this context is a function on T^E. It is not a differential form: 
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rather, it defines the variational problem 
5j(j k 1 )*(L)dt = 

where 7 is a curve in E, j fc 7 is its prolongation to a curve in T*E, and the 
integral is taken with respect to the standard volume form dt on R. The 
Lagrangian is homogeneous if A X (L) = L, A P (L) = for 2 < p < k. 

Each Lagrangian gives rise to its Hilbert form 

defined on T^ h ~ l E; this 1-form has the same extremals as L [3J and is a 
generalisation of the Hilbert form used in Finsler geometry, where k = 1 and 
the formula simplifies to di = SdL. If L is a homogeneous Lagrangian then 
$l is projectable to J 2fc_1 (£', 1). If L (perhaps defined on an open submanifold 
of T k E) has been obtained from a Lagrangian form on some J k n by means of 
the homogenisation trick then the projected image of is just the Poincare- 
Cartan form of the original Lagrangian. 

Finally, the Lagrangian and its Hilbert form may be used together to construct 
the Euler-Lagrange form el using the formula 

e l = dL - d T $L ; 
we therefore see, explicitly, that 

k (_-[\p 

e L = Y, K -^d^dL. 



3 The homogeneous variational bicomplex 

The cohomological approach to the study of variational problems involves con- 
structing suitable sequences of differential forms, or alternatively a bicomplex 
of differential forms. Historically the bicomplex appeared first, and involved 
forms defined on an infinite jet manifold [ 1 , 8 . 9fTU|TT] . The use of a single se- 
quence appeared a little later, and used forms defined on a finite-order jet 
manifold [6]. Subsequently, versions of the bicomplex on finite-order jet mani- 
folds were also obtained [12]. All these constructs were defined initially for jets 
of local sections of a fibration; versions involving contact elements have also 
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appeared, and in some cases the elements of the spaces have been equivalence 
classes of forms rather than individual forms. 



The version of the bicomplex suitable for homogeneous problems was described 
in [TJ, and in its full generality uses certain spaces of vector- valued forms 
rather than scalar forms. But for single-integral variational problems the target 
vector spaces are all one-dimensional, and so we may consider the vector- 
valued forms as identified with scalar forms; that is the approach adopted here. 
The manifolds on which the forms are defined £1X6 ? clS before, the manifolds of 
non-degenerate higher-order tangent vectors T^E on some given manifold E 
with dimE — n+ 1; we write Q r k for the space of r-forms on T^E. (Later on 
we shall talk about certain operators being locally exact, and then it might 
be more appropriate to regard Q r k as the sheaf of germs of r-forms; we leave 
the reader to make this interpretation as required.) 

The basic idea is to combine the usual de Rham sequence of forms with the 
total time derivative operator dy, as in Figure 1. (It is convenient here to 
modify the start of the sequence, replacing the usual 0^M^fi°by0^f2 , 
where O = ft°/R.) In this dia gram we have been explicit about the order 
of the forms; but in fact the operator used to demonstrate exactness of the 
columns will, in general, increase the order: indeed, 'exactness' there must be 
understood as being given modulo a pull-back to a higher-order manifold, and 
so we should imagine not a single bicomplex but a whole family, indexed by 
order. We shall therefore tend to omit any reference to order in the main body 
of the paper, and so improve the clarity of the formulae. In any particular case, 
one may readily calculate the order (or, more accurately, an upper bound on 
the order) of the forms concerned; the general principle is to 'keep going until 
you run out of dots'. 

Our first task is to investigate the exactness of this bicomplex. Define the 
differential operator P : Q r — > Q r (with r > 1) by 



where the sum ranges from zero to a suitably large number so that S p+1 8 
vanishes: if one wished to be explicit one would take a sum up to r k — 1. 

Theorem 1 For every r-form 9 e Q r 
Pd^O = 9 modulo a pull-back. 

PROOF. We use the commutation relation 
Sd^O = d-^SO + r6 modulo a pull-back 



(-1) 



p 



d p T S p 
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Fig. 1. The homogeneous variational bicomplex 

which may be established easily in coordinates for a 1-form and then extended 
to an r-form by induction. Thus 

s p+1 d T e = d T s p+1 e + r{ P + i)s p e , 



and so 



d p T S p+1 d T 6 = (P T +1 S p+1 e + rip + l)d p T S p 9 . 



The coefficients in the alternating sum are chosen so that successive terms will 
cancel out, leaving just the second part of the term for p = which gives 6 as 
required. □ 



It follows from this that the columns of the bicomplex, apart perhaps from 
the first, are globally exact (modulo pull-backs, of course). The first two rows 
are locally exact as they are the usual de Rham sequences; from this we may 
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show by chasing round diagrams that the first column and the bottom row 
are also locally exact. 

Now the bottom row of the bicomplex contains equivalence classes of forms, 
differing by total time derivatives. It is convenient to be able to identify a 
representative form in each class, and the differential operator P allows us to 
do this, at the expense of a possible increase in the order of the forms. 

Lemma 2 For any r-form 9 £ Q r the r-form 

9 - d T P9 

is a canonical representative of the class [9] £ Q r /d,TQ r . 

PROOF. It is clear that 9 — d^P9 £ [9], so we just have to show that start- 
ing with a different form 9 £ [9] we would achieve the same result. But by 
definition 

9 = 9 + d T (j) 

for some <p £ Q r , and 

(dr<f>) - d T P(d T 0) = 

as Pd T (f> = 0. □ 

It is now possible to give a very simple description of the canonical represen- 
tative. 

Theorem 3 Let \& r C Q r (with r > 1) be the subspace of r -forms defined by 
W = {9 £ Q r : S9 = 0} . 

Then id— d^P is a projection operator Q r — > ^> r , so that 9 £ Q r is the canonical 
representative of its class [9] if, and only if, 9 £ \l/ r . 

PROOF. It is an immediate consequence of Lemma [2] that id — d^P is a 
projection operator on Q r , so we need to confirm that its image is \l/ r . From 
the formula for P we see that if 9 £ ^ r then P9 = 0, so that 9 = (id - d T P)9 
and therefore that \P r C Im(id — d^P). 
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For the converse, suppose that 9 = (id — d^P)^ for some G Q r 1 . We note 
first that 



Sd p T = d p T +1 S + r{p+l)d p T 



from the commutation relation mentioned earlier, so that 



Sd T P(^ = Sd T Y, 

v 




d p T S p+1 <j) 




(d p T +1 S p+2 + r(p + l)d p T S p+1 )<j) 



= S<f> 



as the alternating sum collapses once again. We conclude that 
SO = S(<f> - d T P(f)) = S<f)-S<f) = 

so that 9 G \I/ r as required. □ 

We note that, in the case r = 1, S9 = exactly when 9 is horizontal over the 
manifold E. This is not, of course, true when r > 2. 

4 The variational derivative 

The properties of the operator id — d^P suggest that we should consider 
replacing the de Rham differential d : Q r — > Q r+1 by a new operator 

5 : Q r -> fi r+1 , 59 = d9- d T Pd9 

which we shall call the variational derivative. The new operator <5 is clearly 
R-linear, although — despite its name — it is not in fact a derivation. It is, 
however, a coboundary operator. 

Lemma 4 The operator 5 satisfies 5 2 = 0. 
PROOF. From d 2 = and Pd T = id we obtain 
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(d - d T Pd)(d - d T Pd) = -dd T Pd + d T Pdd T Pd 

= -d T dPd + d T Pd T dPd 
= 0. 

The operator 8 mapping Q r to Q r+1 is not locally exact for any r, even modulo 
pullbacks. But by construction 5 takes its values in W +1 C Q r+1 , so we may 
consider the restriction of 5 to \I/ r : 

H° -i» 1 -i» ^ 2 -i» ^ ^ r -i» ^ r+1 -2> . 

This new sequence is indeed locally exact. To prove this, we use the homotopy 
formula for the de Rham differential 

hd6 + dh6 = 6 
where h is the usual Poincare operator. 

Theorem 5 // the r-form 9 E ^ r (with r > 1 ) satisfies 59 = then there is 
an (r — l)-form 



satisfying 5ip = 9. 

PROOF. From 59 = we have 
d9 = d T Pd9 , 

so that 

d T dPd9 = dd T Pd9 = d 2 9 = 

and therefore 

dPd9 = Pd T dPd9 = . 

So Pd9 is d-closed, and locally we may put k = hPd9 so that 
dn = Pd9 . 




V' 1 (r > 2) 
If (r = 1) 
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We now have 

d6 = d^dn = dd^K , 

so that 9 — d^K is (i-closed and locally we may put = h{6 — d^n) G fT -1 so 
that 

d(p — 9 — d^n . 

Then 

S<f> — d(f> — d^Pdcf) 

= (9- d T K) - d T P(9 - d T K) 
= 9 - d T K - d T P0 + d T Pd T K 
= 9- d T Pd 
= 9 

where the final equality arises because 9 G \l/ r is its own canonical representa- 
tive. 

We now consider separately the cases r = 1 and r > 2. If 9 G SI 1 , so that 
G we simply take ip — [(/)] G and then immediately <5-0 = 9. If instead 
r > 2 then G il r_1 , so we put -0 = — d^P<j) G ^ r_1 , and we see that 

5ijj = 5(0 - rf T P0) 

= (id - d T P)d((j) - d T P(f)) 

= d<p — d^Pdcj) — dd^Pip + d^Pdd^Pfy 

= d(p- d T Pd(fi - d T dP(fi + d T Pd T dP(f> 

= dt/>- d T Pd<f) 

= 9 

as required. □ 



5 What has all this to do with the calculus of variations? 

The last two sections have been quite abstract, and we need to return to 
more concrete questions. We do this by relating the first three columns of 
the homogeneous variational bicomplex to specific aspects of the calculus of 
variations. 
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So consider a variational problem S(J Ldt) = 0, for some Lagrangian function 
L G T fc . (The Lagrangian might indeed be homogeneous, but for the purposes 
of the present discussion we do not require this.) We then have [L\ G , and 
we may consider the variational derivative operator defined above. We have 

5[L] = dL- d T PdL 

= dL-d T y ^~ l ^„ d p ^S v+l dL 
p (P + l)! 

= y tl)L d p'sp'dL 

where we have relabelled the sum with p' — p + 1; thus 8[L] = el, the Euler- 
Lagrange form for L described in Section [2j It immediate from Theorem [3] 
that Sel = 0, so that e l is horizontal over E. 

We can also use the local exactness of S to go backwards: given a 1-form ££$', 
when does this represent the Euler-Lagrange form of some Lagrangian? This 
is the 'simple' version of the Inverse Problem of the Calculus of Variations, 
where the multiplier matrix is fixed, and (as in the fibred case) there is a 
definite answer. We consider the Helmholtz-Sonin form 5e G \1/ 2 , and if this 
vanishes then we will be able to construct a (local) Lagrangian for e using the 
formula from Theorem [51 which we may write explicitly as 

L = h(s — dyfiPdE) . 

The question of whether this Lagrangian has minimal order will depend on 
the choice of form k satisfying dn = Pds. 

As an example, we consider the particular case of a second-order form e G 
where in coordinates we have e = Eidq % \ we may think of the form e as a 
second-order differential operator taking its values in the vector bundle T*E, 
and the differential equation itself as the kernel of the operator, a submanifold 
of T^E. The coordinate formula for the Helmholtz-Sonin form is then 

^ 1 f d£i ^ Oej ^ dEi ^ d£j \ ^ .j ^ ^ i 
2 \dq j dq l dip dif ) 

and so the vanishing of these coefficients is a necessary and sufficient condition 
for the form to be locally variational. If this is the case we can use the procedure 
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from the previous section to construct a suitable Lagrangian; we finish by 
demonstrating that making a suitable choice when integrating the two-form 
Pde allows us to find a Lagrangian which is first-order. 

We start by using the condition 8e = to tell us about the structure of the 
functions £j. From the term in dq 3 A dq 1 we have 

dei dej 
dq 3 <9<f ' 



and then from the term in dq 3 A dq l we have just a single term involving q k 3 ^, 

k dej 
q(3) dq k dq 3 



so we see that Si is affine in the second derivative coordinates. Put 
Si = Aij'q J + Bi 



where and Bi depend only on q k , q k and where, from above, = Aji. 
We now consider the construction of a Lagrangian. First, we see that 



de ■ de ■ 

Pde = \^dq 3 A dq l + —^dq J A dq l 
2 dqi dq 3 

i (zk dAjk dBi 

dq 3 dq 3 



and we know that this second-order 2-form must be enclosed. Put Pde = dn, 
where in general 

k = fidq 1 + gidq 1 + 



will be a second-order 1-form, so that 



dn = ^dq 3 A dq 1 + ^-dq> A dq { + ^fidq j A dq { 
aqi oqi ' dq 3 

+ ^dq> A def + '^-dqi A <tf + '^-dq 3 A dg* 
ag J dq 3 dq 3 

+ ^dq 3 A + ^dq> A + ^dq 3 A d<f 
ag J ag- 7 oq 3 
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Of course we are free to add to k any exact 1-form d/j,, where // is a function. 
But from dn = Pde we have 



dhj _ dhj_ _ 
dqi dq l ~ ' 



so we see that there is a function h such that 



— 



We may now consider k — k — dh, and again Pde = dk. We have 



(dq dh \ ( dh dh dh 

fidq 1 + T^dq 1 + — -do 1 - —dq 1 + —-dq 1 + — -dq 1 
dq % dq 1 J \dq l dq 1 dq 1 

\ dq l J \dq l dq 1 J 

= fidq* + 

so that 



dk = ^dq j A dq' + ^dq j A dq' + ^dq> A dq* 
dqi dqi dq> 

+^i dq j A d< ji + ^l d qi A d $ + A dq i . 

ag J ag^ ag-? 

Using the structure of Pde again, we see that /j and (?j are first-order functions. 
Furthermore, from 

®ML - d -h = o 



we see that there is a function g such that 



* ~ dq* 



and we may now consider k = k — dg with, once again, Pde = dk. We now 
have 
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-(*-S)* 

= hdq l 
so that 

d£ = ?J±dq j A rfg 1 + ^dq> A dg l 
ag-? ' aq 3 



Thus 

df 



- — A-- 

dqi ~ 13 



and as = A^ it follows that there is a first-order function / such that 

h dqi 
and therefore that 

10 dq 1 dq j 

We now have a suitable 1-form k, and we know from the theory that e — d^k 
is d-closed, so we can find a function L such that dL = e — d T k. Now 



d T k = d T dq l + ^-dq l 

\ aq 1 J aq 1 



df | ■ df \ dqi+ dl d .i 
dqi dq 1 dq j dq 1 ) dq { 



so that 



, / .j df ., df \ j Of 
e - d T K = \£i- q \ - q \.. — dq - —dq 



dqi dq 1 dqi dq 1 



dq 1 



dqi dq 1 
dq 1 dq 1 



dqi dq 1 
dq 1 - T^-dq 1 



<9<f 
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Thus we see that L is a first-order function, and that 
f = F-L 

for some function F depending only on the coordinates q l . It follows that 

df dL 

dqi dcf dqi dq 1 ' 

and so finally we obtain 

dL 
dcf dqi 
8L_ _ ■ dL 
dq i ^ dqi dq 1 ' 
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